Quantum dynamics and transport in a double well system 
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The simplest one-dimensional model for the studying of non-trivial geometrical effects is a ring 
shaped device which is formed by joining two arms. We explore the possibility to model such a 
system as a two level system (TLS). Of particular interest is the analysis of quantum stirring, where 
it is not evident that the topology is properly reflected within the framework of the TLS modeling. 
On the technical side we provide a practical "neighboring level" approximation for the analysis of 
such quantum devices, which remains valid even if the TLS modeling does not apply. 



I. INTRODUCTION 

In this paper we explore the possibility to model a 
ring shaped device (Fig. 1(a)) as a two level system (TLS) 
(Fig. 1(b)). We shall see that both technical and concep- 
tual difficulties are involved. The model Hamiltonian is 



n = — p2 + Va{x - xa) + VBix- 
zm 



XB 



(1) 



with periodic boundary conditions over x € [— L/2, L/2] 
so as to have a ring geometry, as illustrated in Fig. 1(a). 
Va and Vb represent high barriers, such that the ring is 
composed of two weakly coupled arms. We assume that 
we have control over some geometrical parameters of the 
model and in particular over the heights Xa and Xb of 
both barriers. Our main interest is in the current that 
flows in the system. The current through an arbitrary 
point Xq is obtained as the expectation value of the op- 
erator 

I ^ -—{p5{x-xo) + 5{x-xo)p) (2) 

zm 

Having defined the system and its observables we can 
consider various dynamical scenarios such as coherent 
Block oscillations between the two arms. Then we can 
ask whether a TLS modeling is meaningful. Of particu- 
lar interest for us is the analysis of quantum stirring [l|: 
this means to induce a circulating current by periodic 
modulation of the potential. 

We note that transport due to periodic modulations of 
the potential [2] has been studied mainly in the context 
of quantum pumping [SiiiiSj j where the current is induced 
between reservoirs. The notion of quantum stirring re- 
lates to closed geometry, where the emerging physical 
picture is significantly different (gQ. 

The quantum stirring problem highlights an obvious 
topological subtlety: one wonders whether the non trivial 
topology of the ring is properly reflected in the effective 
TLS model. 

On the technical side we define the unperturbed Hamil- 
tonian Tin as the Xa = Xb = oo limit. In this limit the 
two arms are disconnected from each other, and the diag- 
onalization gives a set of eigen-energies Ei such that each 
eigenstate belongs to only one of the two arms. Then we 
make either Xa or Xb or both finite, and we ask what is 



the perturbation matrix Wij in the reduced Hamiltonian 



TCij = 



El 
E2 



W21 



W12 
W22 



(3) 



Obviously we would like to express the perturbation us- 
ing the transmission coefficients of the barriers. 



II. OUTLINE 



In the ffist part of the paper (Sections IIIHlV|) we es- 
tablish the building blocks. We derive expressions for 
the perturbation matrix Wij and for the reduced current 
operator X^ , and figure out how the topology is reflected 
in the reduced description. 

In the second part of the paper (Sections IVHIVIIj) we 
turn to the applications. We begin with the simplest 
problems: The coherent Bloch oscillations of a particle 
in a mirror symmetric device, and the Wigner decay of a 
particle from a short arm to a long arm. Then we con- 
tinue with the quantum stirring problem, and show how 
one can derive expressions for the geometric conductance. 

In the third part of the paper (Sections IVIIHlX| we 
address some non-trivial technical points that are asso- 
ciated with the analysis, thus exploring the limitations 
of the TLS modeling. We demonstrate that even if the 
TLS modeling does not apply, we still can use a neigh- 
boring level approximation in order to extract results for 
the geometric conductance. 

In the Summary (Section IXIII|1 we highlight the prac- 
tical value of our findings for the purpose of design and 
analysis of quantum stirring devices, and we briefly relate 
to the experimental measurement issue. 



III. THE TLS MODELING SCHEME 

The unperturbed eigenstates tp'^{x) are labeled as 
i — 1,2, corresponding to the two arms of the ring. The 
associated eigen-energies are Ei — kf /{2m). We have 



V'(^)(a;) 



-^ sin(fcia: + (fi) if a; G 1st arm ,.-. 
if a; G 2nd arm 



and a similar expression for ip^'^'{x), where Li is the 
length of the zth arm, and \ip\ < tt/2. Two rep- 



resentative eigenstates are illustrated in Fig. 2. Note 
that the wavenumber of the particle in the ith arm 
is ki = (jr/Li) X integer. Our interest is in a very 
small energy range Ei ^^ E2 ^ E, where the wavenum- 
bers are fci ~ ^2 ~ fc^;, corresponding to the velocity 
VE = {2E/my/^. We would like to ignore all the other 
levels. Later we discuss the validity conditions for this 
TLS modeling scheme. 

Once we lower from infinity one barrier, say barrier A, 
the two states become coupled. In section|IV]we consider 
a delta barrier and obtain the following expression for the 
perturbation matrix: 



ve 



2\/ i^ii^i 



(5) 



where gA is the transmission of the barrier. For i ^ j 
the minus sign is a convention that fixes the gauge (see 
Appendix A). In section |V] we show that essentially the 
same result applies to any other type of barrier, but the 
i = i expression for the energy shift should be somewhat 
generalized. 

If both barriers are finite the two associated perturba- 
tion terms should be added together, and one obtains for 
the energy difference 

£ = (^1 + W^ij + W^n) - {E2 + W^2 + ^§2) (6) 



and for the coupling 



Ve 



2 12 12 2vrrL^ 



(VffI±^V5^) (7) 



The latter expression involves a relative sign zlf that can- 
not be gauged away (see Appendix A). If we had mag- 
netic flux penetrating through the ring we could have, 
instead of the it?, an arbitrary phase factor. 

Using the Pauli matrices we can write the TLS Hamil- 
tonian as 



H, 



n 



(8) 



Defining 6 as the angle between il and the "z" axis, with 
the convention < < tt, the eigenstates hq and ttiq of 
this Hamiltonian are: 



Fo) 



T sin (0/2) \ 
cos(6'/2) J 



m-o) 



cos (0/2) \ 
± sin (0/2) j ^^> 



where the ± indicates the sign of k. The energy difference 
between these eigenstates is 



n = ^/'e 



(10) 



If we have a symmetric well then the effective coupling 
between odd and even levels vanishes (k = 0), and then 
we can get a degeneracy provided we tune appropriately 
the energy level difference e. 

The TLS description is valid if W12 is much smaller 
compared with the level spacing, namely 



max{gA, gs} < L2/L1 



(11) 



where without loss of generality we assume Li > £2- 

In section IIXI we are going to derive expressions for 
the current X^ through barrier A, as defined by Eq.Q 
with xq = xa ■ One observes that the matrix elements of 
this operator in the "standard basis" of Eq.Q vanish, 
because the unperturbed wavcfunctions are zero at the 
barriers. The more careful treatment reveals that the 
reduced operator that gives the net current from the first 
arm to the second arm is 



2" 



V 



(12) 



and it turns out that Xf-: — XaIiji and I^, = Xs^ij, 
where the splitting ratio is defined as 



Xa 



w.i 



W,i + W,^ 



gA^ Jgs 



(13) 



with a similar definition for A^. We have A^ + A^ = 1, 
but contrary to the naive point of view < A^ < 1 is not 
implied. Rather, if the two states have opposite parity, 
then one A is larger than 100%, while the other A is nega- 
tive. We shall see later in Section IXIIII that the physical 
interpretation of the "splitting ratio" requires recogni- 
tion in the existence of induced circulating current in the 
system. Thus the multiple path topology of the system 
is reflected in the TLS modeling via A. 



IV. THE EXPRESSION FOR W^J FOR A DELTA 
BARRIER 

Let us assume that barrier B is infinitely high, while 
barrier A is modeled as a delta function. In other words: 
we consider the simplest possibility of having an infinite 
well [(— i/2) < X < (i/2)] which is divided at x — xa by 
a delta function: 



Va{x - Xa) = XaS{x - Xa) 



(14) 



The total perturbation is obtained from a sequence of 
infinitesimal variations of the barrier height 



n{XA) = ^M-y^ ^dX (15) 

= nioo) + W^ (16) 



For any value of X the Hilbert space of the system is 
spanned by a set of (real) eigenfunction labeled by n. 
The matrix elements for an infinitesimal variation of the 
barrier height is 



dn 

dX 



^'^"^(xa) ^^'"^(xa) 



{m)i 



(17) 



Using the matching conditions for a delta potential at 
X = XA we can express the wave function by its deriva- 
tive: 



r^^'^HxA) 



1 



2mXA 



aV'^"^ (xA+O) - SV'^"^ (xA-O) (18) 



in) I 



A more elegant way of writing this relation is 
1 



V''"n^A) 



2mXA 



J2 dai^^"\xA) 



(19) 



where da is defined as the radial derivative in the di- 
rection of the ath arms that stretch out of the junc- 
tion at cc = xa- Defining the total radial derivative as 
d = di + 82 we get 



dn 

dX 



1 



d-ip'-''\xA) d^'^"'\xA) (20) 



. (2mX^)2 
For a large barrier with small transmission 



9A 



f VE_ 

\Xa 



< 1 



(21) 



the nth and rrith states remain similar to some unper- 
turbed ith and jth states. Accordingly, upon integration 
we get from Eq. lTSl) the result 



Wl 



1 



4m2XA 



dij^^'^XA)] \d^lj^^\xA 



(22) 



Note that in the last equation the contribution to the 
total derivative d comes from one term only, because each 
unperturbed wavefunction t/)'*' [x) is non-zero only in one 
box. Using Eq.dni) we get Eq.©. 



V. THE EXPRESSION FOR Wi^ FOR A 
GENERAL BARRIER 

It is possible to deduce an expression for Wij without 
assuming a specific form of potential barrier. For the pur- 
pose of this calculation we describe the barrier at x = xa 
by a general scattering matrix 



S = e^' 






-V9 



VI -5 c" 



(23) 



Regarding the barrier as a junction it can be embedded 
either in a closed ring geometry with the two arms at- 
tached, or in an open one-dimensional geometry with two 
infinite leads attached. In both cases the differential rep- 
resentation of W should be the same, because W is local 
in space. In other words W^ should come out the same 
for the wavefunctions ^'■^\x) and V'^"'''(^) of the ring, if 
in the vicinity of x = xa they are identical with ^'''(x) 
and ^(^^(x) of the scattering geometry. 

In the scattering geometry it is conventional to label 
the two leads by a = 1,2 and to define a radial coordinate 
r — \x ~ xa\- The flux normalized scattering states of 
the junction (assuming outgoing waves) are \I>(*+). By 
definition we have 



* 



(H 



^—ik^^ 



1 
/ve 



Sue 
[-S2ie'''^'- 



ikET 



] if r G 1st lead 
if r e 2nd lead 



(24) 



A similar expression holds for \I>(2+). If the leads are not 
coupled, the scattering matrix becomes Sq with g = Q. 
In the vicinity of x — xa the unperturbed scattering 
states coincide with those of Eq. ^ up to normalization. 
Namely, in the vicinity of x = xa we have the relation 



vi/W(^) = _^f^ 

VE 



1/2 



pi^i ,l,{i) 



il}'^'>{x) 



where 



'ft 



i(7o + f±ao) 



(25) 



(26) 



with ± sign for i — 1,2 respectively. 

The relation between the scattering matrix and the 
perturbation matrix W can be deduced via the T ma- 
trix formalism. The S matrix is related to the T matrix 
through S — (l — iT)So, or more explicitly 

[SS^%^S,,-^{¥■^^\T\¥^^) (27) 

In leading order T equals W so we have 

($W|T4^|$0")) « i{S^So)So^ (28) 

where 



S-So= e'^° 
Thus 



V~9 



V9 

e-'°"> {S"f - 5a) 



(29) 



Using Eg. ips)) we deduce that each element of 
{tp^'^\W\ip'^^^) involves multiplication by vs/i'iULjy^^, 
while the ao is canceled out. This leads to Eq.® for 
the i ^ j coupling, and a generalized expression for the 
energy level shifts. 



VI. WIGNER DECAY AND BLOCK 
OSCILLATIONS 

If the two arms of the ring have exactly the same length 
Li = L2 = L/2, then the coherent Bloch oscillations of a 
wavepacket in such a symmetric double well are charac- 
terized by the frequency 



n^ 



2\Wi2 



VE_ 

Li 



'9b\ 



(31) 



If one arm of the ring (Li) is very long, and the other 
arm (L2) is short, then a particle placed initially at the 
short arm will decay into the quasi continuum of the long 
arm. The decay rate is given by the Fermi golden rule 



-I 
A 



W12? 



(32) 



where A — {'k/L)ve is the mean level spacing. If the 
arms are coupled through barrier A, while barrier B is 
infinitely high, then the decay rate is 



VE_ 
2L2 



9A 



(33) 



This result agrees with the well known Gamow formula: 
the decay rate is given by the attempt frequency multi- 
plied by the probability to cross the barrier. 

If both barriers are finite, it is important to notice 
that the quasi-continuum of the long arm is composed of 
odd and even states. The state of the short arm, which is 
either even or odd, is coupled to states of the same parity 
with a plus sign in the expression of Eq. ([7]) and to states 
of the opposite parity with a minus sign. Accordingly, the 
decay rate is the sum of the decay rate to states of the 
same parity and the decay rate to states of the opposite 
parity 



Y^ 



27r 



VE 



2L-2 



VE 



2\/LiL2 



9 A ^ y/gs) 



{9 A + 98) 



(34) 



where A± = 2A is the mean level spacing for states with 
the same parity. So inspite of the parity considerations 
we still get the naive result that agree with Gamow for- 
mula. 



VII. QUANTUM STIRRING 

We assume that we have control over geometrical pa- 
rameters of the device, such as the potential floor in each 
arm, the barriers heights, their location, or any other 
gate controlled feature of the potential landscape. With 
a control parameter X we associate a generalized force 
operator 



T 



m 

dX 



(35) 



Quantum stirring means to induce a circulating current 
by changing the parameter X. We assume that the para- 
metric variation is adiabatic so we have a linear relation 
(/) = —GX, where G is know as the geometric conduc- 
tance Q. The Kubo formalism implies that G equals to 
the Berry curvature |9| [l(]||ll| : 



G- E 



m{^n) 



2 lm[Inm]^n 

{Em — En)' 



(36) 



where n is the level in which the particle is prepared, and 
TO are the other levels. 

Within the framework of the TLS modeling the sum in 
Eq. pop contains only one term which involves the states 
no and ttlq of Eq.([9]). For the matrix element of the 
current operator we get 



A— cr„ 



i\- 



(37) 



where A is the appropriate splitting ratio. The matrix 
element of the generalized force operator is calculated 
using Eqs.® and (l35|) 



J~rr 



ds 

dx' 



Ok 
dX 



or. 



-\ 


sin(6l) 


de 
dX 


1 


^ de 




2Vl 


,"dx 


- f 



- cosiO) — — 
2 ^ ' dX 



dX 



(38) 

(39) 
(40) 



where we used ± sin(6') = n/fl and cos(6') = e/fl. This 
leads for the following result for the geometric conduc- 
tance: 



G = 



Xk 
2fF 



de 



Ok 



"dx ax 



(41) 



In the analysis of the operation of a stirring device we 
typically have a well defined region where the potential 
is being varied. We may call this segment "the pump" . 
It is convenient to measure the current elsewhere, where 
the potential is fixed. If barrier A is not part of the 
"pump" then we can measure the current at xo = xa- 
Then it follows from the definitions of A and k that the 
product Xk does not change with time, even if barrier B 
is modulated. Then we can rewrite the above formula as 



G 



XqKq 



de 



Ok 



'dx dx 



(42) 



where Ao and kq are that values at some arbitrary mo- 
ment of time. Typically the variation of X leads to a 
sequence of level crossings if k is disregarded. These be- 
come avoided crossings if k is taken into account. At 
the vicinity of a crossing we typically can use a linear 
approximation: 



e 

K 



ex{X- Xo) 

kq + kx {X - Xq) 



(43) 
(44) 



The amount of probability dQ — Idt which is being trans- 
ported equals —GdX. For an individual crossing the dX 
integration over G can be performed using 



-l-CXD 



a{b + ex) dx 



2a 



-00 (a2a;2 + (& + ca;)2)3/2 bVa^ + c^ 
+°° caxdx 2c2 



-00 (a2x2 + (fe + ca;)2)^/2 abVa^ + c^ 

Then we get the result 



Q 



±x 



o^l + {k/ef 



(45) 



(46) 



(47) 



where the ± is determined according to the sign of e. We 
observed that in order to get the "quantized" value Q = 1 
there should be neither topological splitting (A = 1) nor 
barrier modulation [k = 0) during the transition. 



VIII. THE NEIGHBORING LEVEL 
APPROXIMATION SCHEME 

A major interest is in systems with zero temperature 
Fermi occupation. In such a case Eg. ([5^ has to be 
summed over n up to the Fermi energy. It turns out 
[l| that the result is dominated by the contribution that 
come from the coupling between the last occupied level 
and its neighboring empty level. This suggests to adopt 
a neighboring level approximation scheme that holds ir- 
respective of the validity of the TLS modeling, and coin- 
cides with it if the condition of Eq. (fTT|) is satisfied. The 
key idea is to characterize each eigenstate by a mixing 
parameter 



while the true eigenstates are with (see Fig. 2) 



e 



2 arctan 



/Prob(a; G 2) 
Prob(a; G 1) 



(48) 



such that 8 = for states that belong to the first arm 
and 8 = TT for states that belong to the second arm. 
Numerical examples are presented in Figs. 3-4. If we are 
given 8 then we can construct the eigenstate using a 
procedure that we describe below. If the TLS modeling 
applies then 8 becomes essentially the same as 6. 

Let us see how we construct the wavefunction given 
the energy E = £'„, the mixing parameter 8 = 8„, and 
the parity -^ with respect to (say) barrier A, as defined 
in Appendix A. Consequently it is convenient to set the 
origin such that xa = 0, and write the nth eigenstate of 
the ring as 



V'*"^^;) 



zt?Ci sin(fcia; -I- <y5i) if x e 1st arm 
C2 &va{k2X + (^2) if a: G 2nd arm 



(49) 



where Ci > 0, and \ip\ < 7r/2. Assuming ksL ^ 1, the 
amplitudes satisfy the normalization condition 



h.c,^ 


+ ^L2C2^ « ] 


follows that 




Ci « 


nr /e 

\ / -;— cos — - 


C2 « 


nr . /8 

Vi2"H2 



1 



(50) 



(51) 



(52) 



We still have to say what are the wavenumbers ki and 
k2, and the phase shifts ip. Let us see first how they are 
determined within the framework of the TLS modeling, 
and then how they can be found irrespective of the TLS 
modeling. 

Naively the |no) and |mo) eigenstates, within the 
framework of the TLS modeling, are the superposition 
of the basis states of Eq. ((3]) and accordingly 

8 = 61, -K-e (53) 

fci = corresponds to the unperturbed Ei (54) 

fc2 = corresponds to the unperturbed E2 (55) 

ipi — same as the unperturbed (56) 

ip2 = same as the unperturbed (57) 



8 


w 61, TT-e 


(58) 


ki 


= corresponds to En 


(59) 


k2 


— corresponds to En 


(60) 


^1 


— shifted 


(61) 


^2 


= shifted 


(62) 



To be more specific, we have S^™"-* « 6 and 8^""^ « tt—O 
and hence 8('"«)-H8("o) « tt if the TLS modehng is valid 
(see Fig. 4). We note that from Eq. (|5HI52[) it follows that 
within the framework of the TLS approximation we have 



p("io)^(no) 



— sin(6l) 
Li 



(63) 



This will be used later on in order to obtain simplified 
expressions for the matrix elements of various operators. 

Whether ki and k2 in Egs. (1541155]) correspond to the 
same energy or not, is not a big difference for us because 
we assume /ci ~ fc2 ~ ks in any case. The main problem 
with the naive version is related to the phase shifts, as 
demonstrated in Fig. 2. The variation of the phase shift 
as the barriers are lowered reflect that there is a non-zero 
probability to find the particle in the region of the barri- 
ers. In particular if the phases ipi remained the same it 
would imply that all the matrix elements of J"* and I^ 
would be zero. It is essential to take the variation of (p 
into account in order to get a non-zero result for the geo- 
metrical conductance. We shall discuss the calculation of 
the matrix elements Inm and J-'nm in the next sections. 
First we would like to discuss how the required informa- 
tion on the variation of the phases ip can be extracted. 

In order to express (/s by 8, we write the wave function 
of Eq. (l49| as ingoing and outgoing waves and set the 
origin a; = at either one of the barriers, for example 
barrier A. We match the wave functions of the two bonds 
by the barrier scattering matrix 



C2 e+'^^ 



- S, 



^Ci e- 



lipi 



2 e 



(64) 



and get closed equations for the phase shifts 



yr^sin(2^i-a-7)=l-|(l+(^) ) (65) 



Vl^gsin(2(y92 +a- "f)^l --^11 




(66) 



So far everything is exact. So once we have 8 we 
can find the phases and construct the wavefunction. We 
would like to focus in the rest of this section in the regime 
where the TLS modeling applies. Assuming that 8 is de- 
termined by 9 we want to find what are (pi and (p2, so 
as to construct a proper wavefunction. Neglecting terms 
of order g and expanding arcsin(l— .t) as n/2 ± ^/2x we 
obtain 

7 + a TT v^ [L^ . f'd\ if.^. 



where the ilf sign should be the same as in Eq. (|64| , which 
can be estabhshed by direct substitution. A similar ex- 
pression can be obtained for Lp2- We note that within the 
framework of this TLS approximation we have 



^^""^ 



^("o) 



±\/5l 



1/2 



sin(6') 



(68) 



where the sign is the same as that of k. We now have 
all the building blocks needed for the calculation of the 
matrix elements. 



X. STIRRING BY BARRIER MODULATION 

In this section we calculate the geometric conductance 
as determined by the matrix elements of the generalized 
force that is associated with modulation of a delta bar- 
rier. The motivation is to verify the results of the reduced 
description against the direct full Hilbert space calcula- 
tion. The potential barrier is given by 



Vb{x) = Xb5{x - xb) 



(74) 



IX. THE EXPRESSION FOR J„ 



The stirring is induced by variation of the barrier height 
Xb- The associated generalized force is 



If we adopt the TLS point of view, we can postulate 
a self-consistent definition of the current operator based 
on the continuity equation. For this purpose we define 
the occupation operator A/" for one of the arms as 



N - 



1 




(69) 



and deduce the definition of the current operator from 
d 



dt 



N = i[nM] = ^ 



(70) 



where X is given by Eq. p2|) . If we turn off the coupling 
at barrier A we get the same expression multiplied by 
\b , while if we turn off the coupling at barrier B we get 
the same expression multiplied by A^. 

The above reasoning bypass the confrontation which is 
involved in carrying out a direct calculation, and hence 
contains an uncontrolled error which is associated with 
the assumption that a TLS description of Hilbert space 
is valid. If we revert to the original definition of Eq.©, 
then the matrix elements are given by 

In™ = t:^ (a^<"^ ^("' - V'^"^ a^^"0 (71) 

zm v / x=xo 

For the calculation of I;^„ we set xo = xa = 0. As was 
already pointed out, in order to get a non-trivial result, 
we have to take into account the phase shifts ip which 
was calculated in the previous section. Substituting the 
wave function of Eq. ((49)) we get 



I^ = 

nm 



1 






-i—C]"''C\ 
2m 1 1 






(m) 






■ I (™) 

Sin((y9j^ 



"^1 



(72) 



Whenever the TLS modeling applies we can substitute 
Eqs.dM} and dM]) into Eq.lO. Neglecting the second 
term we get 



1^ 



T« 



Ue 



2^/^^c^ 



(73) 



where the sign is the same as that of — k. One notices 
that the expression for T-^^^^ can be written as Eq. (l571) 
where k and A are given by Eq.([7]) and Eq. (|13p . 



T = - 



dU 



6{x - xb) 



OXb 
with the matrix elements 



(75) 



(76) 



For the wavefunctions amplitudes we use Eq. (j63p and for 
the phase shifts Eq. ([67)) . We also substitute the scatter- 
ing matrix parameters that describe a delta barrier 



IB ~ -7r/2 + ^/g^ 
as = 



(77) 
(78) 



where the approximation is valid for 55 ^ 1 and the 
relation of gB and Xb is given in Eq. (PT|) . With the 
above approximations we get 

^mono - 9B ^'~/'' sin(0) T ^^^ cos(0) (79) 



4L1L2 



2ViIi^ 



where the sign should be the same as that of ^n. In 
order to verify the consistency with the TLS expression, 
we differentiate Eq.® and Eq.©: 



de 

dX^ 
dn 

dX^ 



= 9B- 



L2 — Li 
2L1L2 

9B 



VL1L2 



(80) 
(81) 



and substitute into Eq. ([38l) . Indeed we obtain the same 



result for Tmono as above. 



The geometric conductance of Eq. ([36)) involves the 
multiplication of J-'mono with Tnomoj leading to 



G 



1 2 L2-L1 gAgB^gA^/^gs'/^ 

4''^(LiL2)^ ^' 

1 2 L2 + L1 gA gB + gA^'^gs^'^ 

^ ~ '"ETT~rT2 



{LIL2Y 



fi3 



(82) 



The calculation of the transport proceeds as in Sec- 
tion |VlTl 



XI. STIRRING BY BARRIER TRANSLATION 

In complete analogy with the previous section we 
would like to calculate the geometric conductance as de- 
termined by the matrix elements of the generalized force 
which is associated with the translation of the barrier: 



J" 



m_ 

dXB 



Xb5'{x- xb) 



(83) 



One obtains 



r^^C") ^(™) + ai/;(™) ^(«)^ (84) 



where dip is the average derivative on both sides of the 
barrier. We simplify this expression by using Eq. pT 
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Assuming high barriers we get in leading order 



which together with Ea. ([S5)l leads to 



(m)^(n) 



a 



J~rr 



1 2 



L1+L2 
L1L2 



sin(6i) 



(85) 



(86) 



(87) 



In order to compare the above result for J-mano with the 
TLS result of Eql38]we calculate the variation of the po- 
tential floor by taking in Eq.® the energies of infinite 
wells with Li = xb and L2 ~ L — xb- We get 



de 

dXB 

dn 

dxB 



Li + L 



^^ ~^ — 172 (Vm ^ V9b) 



{L1L2) 



(88) 
(89) 



Substitute into Ea. ([38|) indeed leads to the same result 
for J-'mono as abovc. Note that in this case (unlike the 



previous section) the second term in Eg. psp which in- 
volves the variation of k is of higher order in gB and 
therefore should be excluded. 

The geometric conductance of Eq. ([36|) involves the 



(90) 



multiplication of Tmana with Xn^^rao-, leading to 



G 



1 ^ 4 Li+ia 
2^^^ ^ {L^L2f 



gA^JgB 



n^ 



The calculation of the transport proceeds as in Sec- 
tion IVIII One realizes that a translation of the barrier 
is effectively equivalent to the variation of the potential 
floor difference, as long as it does not involve modulation 
of its transmission (which is assumed to be small). 



XII. ERROR ESTIMATES AND LIMITATIONS 

If we vary a parameter X then the energy levels En{X) 
form a "spaghetti" which is characterized by a mean level 



spacing A and possibly by narrow avoided crossings with 
splitting Aq. For the ring system that we are considering 
it follows from the estimate of k that 



Ao 

A 



— - - min{l, ^/bg} 



(91) 



where h — L1/L2 and g = maxj^^j^s}- The condi- 
tion Eg. pTj) for the applicability of the TLS modeling 
ensures Ag <C A. In such circumstances Eg. ([55]) for the 
geometric conductance, which in essence is a sum of the 
type X]^o(^o "l""^)"^' implies that the error that is 
involved in the neighboring level approximation is 



error(G) 
G 



Ao 
A 



^ - b9 «1 



(92) 



Once the TLS modeling fails the error becomes of order 
unity. This sounds bad, but in fact it is not so bad. 
The good news is that the far levels contribute to G a 
correction which is of the same order as the leading term. 
Therefore with the neighboring level approximation we 
can still get a realistic estimate disregarding numerical 
prefactors of order unity. 

Having 6 ^ 1 is very interesting, because then we have 
a non-trivial intermediate regime 1/h <C <? «C 1 where nei- 
ther 1st order perturbation theory with respect to "zero" 
height barriers, nor 1st order perturbation theory with 
respect to "infinite" barriers applies. This is the regime 
where each level of the small arm forms a distinct Wigner 
resonance with the quasi-continuum states of the long 
arm. Obviously the TLS modeling is not applicable in 
this regime, but the neighboring level approximation still 
provides a decent starting point for a calculation. We 
shall explore this Wigner regime in a future work. 

One may also wonder whether the specific results that 
we have obtained for stirring using a delta barrier ap- 
plies also for a thick barrier. On physical grounds it is 
quite obvious that the induced current is determined by 
the scattering matrix of the modulated barrier. Conse- 
quently if the S{E) of the modulated barrier is E in- 
dependent within the energy range of interest, it can be 
regarded as representing a delta function, and the results 
should come out the same. 

Finally one may wonder about the implications of fi- 
nite temperature or non-adiabatic driving. These aspects 
are complementary to the theme of the present paper. 
Namely, as discussed in [1], at finite temperatures the 
statistics of the occupation should be taken into account. 
So we have to average (so to say) over the level that we 
have labeled as ng with an appropriate weight as im- 
plied by the Fermi function. On the other hand the non- 
adiabatic effects require to introduce in the denomina- 
tor of the Kubo formula Eq. ((36)) a term that represents 
the "width" of the Fermi-golden-rule transitions. Then 
the weight of the neighboring level in the sum becomes 
smaller compared with the total weight of the far levels. 
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XIII. SUMMARY 

We have developed a practical procedure for the anal- 
ysis of a one dimensional double well system, which is 
both powerful and illuminating. The procedure assumes 
that we have a way to find the cign-energies En of the de- 
vice, and the mixing ratio 0„ of each of them. Given the 
transmissions of the barriers we further characterize the 
device by the splitting ratio A. With these ingredients in 
hand we can analyze any stirring process and obtain ex- 
plicit expressions for the geometric conductance G. The 
calculation simplifies if the TLS modeling applies, be- 
cause then the mixing ratio can be determined form the 
diagonalization of a 2 x 2 matrix. 

In particular we obtain explicit expressions for G due 
to either barrier translation (generalizing a result that 
has been obtained in [l|), or barrier modulation (gener- 
alizing a result that has been obtained in 7']), and verify 
that they agree with the naive self-consistent TLS calcu- 
lation. We see that whenever the TLS modeling applies 
the proper calculation in the full Hilbert space gives the 
same result as the naive calculation in the TLS Hilbert 
space. 

As a by product of the TLS analysis wc find that the 
pumped "charge" during an avoided crossing is not quan- 
tized (see Eq. (|47| ). not only because of the topological 
splitting effect, but also due to a dynamical effect that 
arises if the barrier is modulated. 

The practical importance of the TLS modeling in con- 
dense matter physics is obvious. On the other hand the 
specific application to the study of quantum stirring de- 
serves a few words regarding the measurement proce- 
dure and the experimental relevance. As explained in 
[l| it should be clear that the measurement of current 



Ullii 



in a closed circuit requires special techniques 
These techniques are typically used in order to probe per- 
sistent currents, which are zero order (conservative) ef- 
fect, while in the present paper we were discussing driven 
currents, which are a first-order (geometric) effect. It is 
of course also possible to measure the dissipative conduc- 
tance (as in |12||V During the measurement the coupling 
to the system should be small. These are so called weak 
measurement conditions. More ambitious would be to 
measure the counting statistics, i.e. also the second mo- 
ment of Q as discussed in J15llla | which is completely 
analogous to the discussion of noise measurements in 
open systems [17iil8i [. Finally it should be pointed out 
that the formalism above, and hence the results, might 
app ly to experiments with superconducting circuits (see 



The unperturbed eigenstates are labeled by i (or option- 
ally by j). In the TLS scheme z = 1, 2. If the coupling 
is non-zero the exact eigenstates are labeled by n (or op- 
tionally by m) . Within the framework of the neighboring 
level approximation scheme we focus on two levels that 
we label as n = ng and m = ttiq. If the TLS modeling 
applies then the states ng and ttiq are regarded as linear 
combinations of i = 1, 2. 

The unperturbed states ? = 1 , 2 are characterized by 
their parities zi^ and z^ respectively. The relative sign zlf 
in Eq. ([7]) equals the product of zt^ and z^ . Inverting the 
arbitrary gauge sign of either tp^^>{x) or tp^'^'{x) would 
multiply the expression in Eq.® by a global minus sign, 
while the relative sign zlf remains unchanged. The gauge 
invariant relative sign is due to the fact that the unper- 
turbed states are either odd or even: we have plus sign if 
both states have the same parity and minus sign if they 
have opposite parity. 

Each exact eigenfunction n, as written as in Ea. (|49p . 
is characterized by what we call the parity d? with re- 
spect to barrier A. Positive parity means that the radial 
derivatives as defined in Eq. p^ have both the same sign. 
Optionally we can define zi^ as the parity with respect to 
barrier B. This parity ilf is not a symmetry related quan- 
tum number, but it is merely required in order to define 
the wavefunction of Eq. (j49p in a unique way given the en- 
ergy and the mixing ratio. If the TLS modeling applies 
then for positive (negative) k the state uq of Eq. ([9]) has 
negative (positive) parity, while the niQ state has posi- 
tive (negative) parity. Within this framework the parity 
zl^ with respect to barrier B is ilf multiplied by df . 

We have verified that the various ± signs through the 
paper are consistent, which is not always evident in a 
superficial look. 
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APPENDIX A: CONVENTIONS AND 
NOTATIONS 



Consider two segments that are connected at points 
that are labeled as xa and xb- In the absence of cou- 
pling each segment is regarded as a one dimensional box. 
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FIG. 1: Panel (a): Illustration of a ring shaped device that 
is divided by the barriers Va and Vb into two arms of length 
L\ and 1/2- The current is measured through the dashed sec- 
tion near barrier A. In the quantum stirring scenario it is 
assumed that there is a gate control over the potential floor 
of each arm, or over the height or the location of barrier B. 
Panel (b): Within the framework of the TLS modeling, the 
reduced Hilbert space contains two levels. The perturba- 
tion Wij is due to having finite rather than infinite barri- 
ers, so it corresponds to the difference Ti — 7i(oo) and not to 
1/ = 7i — 7i(0). See the text for further details. 





FIG. 2: Upper panel: Two nearly degenerate eigenfunctions 
■i/)(a;) of a particle in a ring with arms of length Li = 1 and 
L2 = 2.23. These are the two unperturbed states of Eq.Q. 
Lower panel: The exact eigenfunctions assuming that the bar- 
riers are finite {qa ~ 0.28 and qb ~ 0.06). These do not van- 
ish at the barriers, and therefore cannot be written as a super- 
position of the unperturbed states. Still we explain in the text 
how a decent approximation for the former can be obtained 
using the neighboring levels approximation scheme. 




FIG. 3: We consider a particle of mass m = 1 in a ring of 
length L = 151.43. The position of barrier B is X, so we 
have L\ = X and L2 = L — X . We calculate numerically fc„ 
and 6'"' for two neighboring levels (solid and dashed lines). 
The sum 9('")-(-9(") is plotted as a dash-dotted line. We have 
high barriers with qa ~ 10~^ and qb ~ 10~^. Accordingly we 
expect TLS modeling to be valid: The dotted lines indicate 
the values O — -k 12 (expected crossing point) and = tt (ex- 
pected sum) . For sake of comparison there is a third dotted 
line that indicates the value of O that corresponds to equal 
amplitudes C\ — €2- 



2.9 




2.895 




2.89 




2.885 


■^ ^ ^^—^ 


2.88 




2.875 




2.87 




2.865 





2.875 

2.87- 
2.865 



FIG. 4: The same as the previous figure, but here the TLS 
modeling does not apply. In the left panels one barrier is high 
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and one barrier is low {qb ~ 0.9), while in the 



right panels both barrier are low [qa ^ Qb ^ 0.9). 



